Instability of Walker Propagating Domain Wall in Magnetic Nanowires 
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Stability of the well-known Walker propagating domain wall (DW) solution of the Landau-Lifshitz- 
Gilbert equation is analytically investigated. Surprisingly, the Walker's rigid body propagating DW 
mode is not stable against the spin wave/wavepacket emission. In the low field region only stern spin 
waves are emitted while both stern and bow waves are generated under high fields. In a high enough 
field, but below the Walker breakdown field, the Walker solution could be convective/absolute 
unstable if the transverse magnetic anisotropy is larger than a critical value, corresponding to a 
significant modification of the DW profile and DW propagating speed. 

PACS numbers: 75.60.Jk, 75.30.Ds, 75.60.Ch, 05.45.-a 



Magnetic domain- wall (DW) propagation in nanowires 
has attracted considerable attention in recent years [lHfa] 
because of its fundamental interest and potential appli- 
cations @, y] ■ Field-driven DW dynamics is governed by 
the Landau-Lifshitz-Gilbcrt (LLG) equation which has 
a well-known Walker's exact rigid-body propagation so- 
lution |l| for a one-dimensional (ID) biaxial wire. This 
Walker solution plays a pivotal role [6-8] in our current 
understanding of both current-driven and field-driven 
DW propagation in magnetic nanowires. A genuine so- 
lution of a physical system must be stable against small 
perturbations. Although there is no proof of the sta- 
bility of the Walker solution and there are signs [9|, [l0( 
that this solution may be unstable, at least under cer- 
tain conditions, the validity of the Walker solution for 
a ID wire is always taken as self-evident. Any devia- 
tion in experiments or numerical simulations are assumed 
to be attributed to the quasi-lD nature or other effects 
[7|. On the other hand, applications of spintronics de- 
vices require accurate description of DW motion [Ill4l4 1 . 
Thus, the stability of the Walker propagating DW solu- 
tion becomes vital in our understanding of DW propaga- 
tion along a magnetic wire. 




FIG. 1: (Color online) Illustration of transverse head-to-head 
DW of width A in a nanowire, with easy axis along z and hard 
axis along x. In the absence of external magnetic field (up- 
per), a static DW exists between two domains with m z = ±1. 
Under a field parallel to the easy axis, the Walker propagating 
DW moves towards the energy minimum state (m z = — 1) at 
a speed v while the DW profile is preserved. 



mensionless ID LLG equation [5J, 
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In this paper, by using a recipe that is based on a 
series of recent advances in nonlinear dynamics theory, 
the stability of the Walker's exact DW solution is the- 
oretically analyzed. To our surprise, the solution is not 
stable against spin wave emission. In the low field re- 
gion, only stern spin wave could be observed while both 
stern and bow waves emerge under high field. Severe 
distortion in propagating DW profile and resulted devi- 
ation of DW speed from the Walker formula can occur 
when the transverse magnetic anisotropy is larger than a 
critical value and the external field is sufficient high, but 
below the Walker breakdown field. For a given transverse 
magnetic anisotropy, the solution is transient unstable at 
low field and convective/absolute unstable at high field, 
corresponding to emitted different spin wavepackets. 

To study the stability of Walker's exact propagating 
DW solution under an external field, we consider the di- 



This LLG equation describes the dynamics of magneti- 
zation M of a magnetic nanowire schematically shown 
in Fig. [1] With the easy axis along the wire (z direc- 
tion) and the width and thickness being smaller than the 
exchange interaction length, exchange interaction domi- 
nates the stray field energy caused by magnetic charges 
on the edges; the DW structure tends to be homoge- 
neous in the transverse direction [15J, i.e., behaves ef- 
fectively ID. We are interested in the behavior of a 
head-to-head DW under an external field shown in Fig. 
Q] In Eq. ([TJ, m is the unit direction of the local 
magnetization M = mM s with saturation magnetiza- 
tion M s and a is the phenomenological Gilbert damp- 
ing constant. The effective field (in the units of M s ) 
is h e ff — Kum z z — K±m x x + Ad 2 m/dz 2 + Hz where 
K\\ , K± , and A are respectively the easy axis anisotropy 
coefficient, the hard axis anisotropy coefficient, and the 
exchange coefficient. H is the external magnetic field 
parallel to z. The time unit is (7-M s ) -1 , where 7 is the 
gyromagnetic ratio. Using polar angle 9 and azimuthal 
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FIG. 2: (Color online) Left are the essential spectrum (shad- 
owed regions) for H = (a) and H « 0.02# c (b). The Fred- 
holm borders are A* 2 (&)- Solid border lines correspond to 
spin waves with negative group velocities while the dashed 
border lines are for the spin waves with positive group ve- 
locities. Propagating DW wall emit stern waves in low 
fields (right of (a)), and stern and bow waves in higher field 
(Cl.02.Hc < H < H c ) (right of (b)). The green dots are zero 
group velocity modes. K± = 1 is used. 



angle ip for m as shown in Fig. [TJ this LLG equation has 
a well known Walker propagating DW solution [l| , 



sin2(p w (z, t) 



H 



In tan —9, 
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vt 
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Here H c = aK±/2 is the Walker breakdown field and 
A = (K\\/A + cos 2 ip w K ± /A)- 1 / 2 is the DW width which 
will be used as the length unit (A = 1) in the analysis 
below, v = AH/a is the Walker rigid-body DW speed 
that is linear in the external field and the DW width, and 
inversely proportional to the Gilbert damping constant. 
Solution $Z§ is exact for H < H c . 

To prove the instability of solution © against spin 
wave emission, we follow a recently developed theory 
(Sandstede, 2001 0; Fiedler and Scheel, 2003 0) for 
stability of a general travelling front such as a propagat- 
ing head-to-head DW shown in Fig. [TJ Consider a small 
deviation of the Walker solution, 9 w + 9 and ip w + ip with 
\9\ <C 1 and \<p\ <C 1, the equations satisfied by 9 and 
ip can be readily obtained from Eq. (UJ. In the moving 
frame of the DW velocity v (with coordinate transforma- 
tions of z — >• £ = z — vt and t — > t), the linearized equa- 
tions of 9 and ip in a two-component form of A = (9, ip) T 
(superscript T means transpose) are 



rfA 
~dt 



= L(9 w , Vw ,d/dtd 2 /d 2 0A. 



(3) 



L is an inhomogeneous operator, depending on £ through 
9 W . The energy dissipation due to time variations of 9 



and (p are absent because they are higher order terms [5j . 
It shall be important when the deviations are out of the 
linear regime. The possible solutions of Eq. ([3]) of type 
A = A 1 (£)e A * define spectrum A of L. Similar to energy 
spectrum of a quantum system, A can be continuum and 
discrete. The former is often called essential spectrum 
while the later point spectrum. The spectrum A shall 
determine the stability of the Walker solution. In terms 
of A' = (9, ip, 89/ d£, 8ip/8£,) T , equation ^ becomes a 
four dimensional first-order ordinary differential system 



-A' = TA', 



r = 



o I, 
B C 



(4) 



here I is a 2 x 2 identity matrix, and 2x2 matri- 
ces, B and C, have following matrix elements: B\\ = 
aX/A + (1/2A)(K X - 2K\\ - K±y/l - p 2 )cos[2G(f)J - 
H tanh£A4 where G(£) is the Gudermannian function 
and p = H/H c ; B 12 = (A - K±ptanhg)/(Aco8h(;); 
B 21 = -cosh£(2A + K_ L pta-nh^)/(2A); B 22 = {aX - 
K ±y /T=?)/A', Ch = -va/A; C 12 = -v/iAcosh^; 
C 21 = vcosh^/A; C 22 = -2 - va/A. 

According to the theory of Rcfs. |16l420| . the spectrum 
of L is determined by Eq. (j4j). The essential spectrum is 
bordered by the Fredholm borders (defined below) of (|4j 
with T replaced by its two limits of £ — > ±oo, denoted as 
r ± = lim^iooF. Since T^ are constant 4x4 matrices, 
solutions of Eq. ((4]) with T — T ± are linear combina- 
tions of Aoe" ^ with k being complex numbers. Pure 
plane wave solutions (k = ik) exist only when A sat- 
isfies det(r ± (A) + ik) =0 with k £ (—00,00). Each of 
the two equations has two branches of allowed A labeled 
as Aj = 2 (fc), known as the Fredholm borders [la. Il7l Il9j. 
In another word, Eq. (gj) with T = T+ (r = T - ) has 
pure plane wave solution when A is on Xf 2 (k) (A 1 " 2 (fc)). 
For those A not on Xf 2 (k), each Eq. (g) with T = T ± 
have four k . s whose real parts are nonzero. If one uses 
(n + , n_) to denotes A for n + (n_ ) being the number 
of « s with positive (negative) real part. Then both 
A^ 2 (fc) an d X^ 2 (k) divide A-plane into three parts with 
(ji + ,n^z) = (1,3), (2,2), and (3,1), respectively. Ac- 
cording to theorem 5.35 in Chapter 4 of Kato [20j, the 
essential spectrum of L must be in the regimes with 
nl + n^_ =/= 4. For A on boundaries A^ 2 , the associated 
eigenmode is plane wave (spin wave) while eigenmodes 
for A not on the boundaries are spin wavepackets. 

In order to understand numerical results in Ref. [l0(, 
parameters of yttrium iron garnet (YIG) [14| are as- 
sumed in our analysis with A — 3.84 x 10 -12 J/m, 



7 = 35.1kHz/(A/m), and M s = 
- 0.001 is used and K± is a vary- 



A' t | = 2 x 10 3 J/m 3 

1.94 x 10 5 A/m. a 

ing parameter. Fig. [2] plots the essential spectrum for 

K± = 1. In the absence of an external field, two branches 

of the spectrum of T ± are the same, Xf 2 (k) = A 1 " 2 (fc), 

shown in Fig. [2ja). Since the spectrum encroaches the 

right half plane, unstable plane waves shall exist and spin 
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FIG. 3: (Color online) (a) A" 2 for K x = 1, p = 0.35 (dashed 
curve) and 0.36 (solid curve). The absolute spectrum is be- 
tween two branching points Sdi and Sd,2 (green dots). Inset: 
Plot of 7?e(fi;j) and Re(ng) vs. A between Sdi and Sd,2- At 
Sdi t 2, ^2 — K 3 ■ 3 ) Graphical illustrations of three types 
of instabilities caused by unstable absolute spectrum. Green 
curves indicate initial profiles of unstable modes while the 
dotted (blue) and dashed (red) curves are their later profiles. 
A transient instability (i) emits unidirectional waves (propa- 
gating to the left). A convective instability (ii) emits waves 
in both directions. An absolute instability (iii) emits waves 
that do not travel in the moving frame, or move with the DW. 
(c) Phase diagram of transient (TI) and absolute/convective 
(AI/CI) instabilities. The boundary is the bifurcation line 
between Tl/CI-and-CI instabilities in Kx and p = H/ H c 
plane. The bifurcation line is only plotted for K± > K± 
here K ± « 0.085 at which H 2 - H c (p = 1). Noted that our 
analysis is valid for fields below the Walker breakdown value. 



wave emission are expected. Similar conclusion was also 
obtained in early study [21|. Solid lines are for negative 
group velocity (determined by Im(dX/dk)), thus these 
are stern modes. The dashed lines indicate positive group 
velocity, corresponding to bow modes. The green dots 
are zero group velocity points. According to Fig. Ufa), 
all unstable modes have negative group velocities so that 
DW can only emit stern waves in the low fields. As the 
external field increases, A+ 2 (fc) and Aj~ 2 (fc) are different, 
and the area of essential spectrum in A-plane becomes 
bigger and bigger (shadowed regimes in Fig. |2Jb). The 
green dots also moves toward 7m(A)-axis and cross it at 
H « 0.02H C (Fig. m»). Further increase of H, the 
unstable modes have both positive and negative group 
velocities although the most of them have the negative 
ones. One shall have propagating DW to emit both stern 



and bow waves. The stern waves should be stronger than 
the bow waves as schematically shown in the right figure 
of Fig. HJb) . This is exactly what were observed in nu- 
merical simulations on dissipative wires for stern wave 
emission in low field [9( and stern-and-bow wave emis- 
sion in high field [lOj. In a realistic wire with damping, 
emitted spin waves will be dissipated after a short dis- 
tance, and are hard to be observed in experiments. 

Interestingly, severe instabilities of a propagating DW 
are not determined by the essential spectrum, but by the 
absolute spectrum |16h201 I22H24J that could change the 



DW profile so that DW propagating speed would be sub- 
stantially modified |5j . As mentioned early, for each A in 
the complex plane, there are four n i (i = 1,2,3,4) for 
T ± , ordered by their real parts as Re(n 1 ) > Re(n 2 ) — 
Re(n 3 ) > Re(n 4 ). Then A is said to belong to the ab- 
solute spectrum if and only if Re(iit(X)) — Re(K 3 (X)) 



or Re(n 2 (A)) = Re(n 3 (A)) 



The branching 



points are special points in the absolute spectrum, de- 



(A sd ) 



z (X sd ). They are 



noted as X s d, satisfying k^ 

non-traveling modes [24j, [27[ . For Kj_ -- 1 , the absolute 
spectrum in the right half A-plane is generated by T~ . 
Fig. [3Ja) shows two branches Aj" 2 . They are well sep- 
arated by the real axis for p = 0.35 as shown in Fig. 
[3ja) (dashed curves) and no absolute spectrum could be 
found in the right half plane. As the field increases, the 
two branches get closer with each other and at an onset 
field H2, depending on Kx, two branches tangent at the 
real axis and then separate again in horizontal direction 
as shown in Fig. EJa) for p = 0.36 (solid curves). At 
this moment, the absolute spectrum begin to emerge on 
the real axis (the segment between two branching points 
5^1,2 (green solid dots). The dependence of ^(k^) or 
Re(K 3 ) on Re(X) between these two points is shown in 
the inset of Fig. [3Ja) (solid segment). 

According to Refs. |22h24| . wavepackets would be 
emitted if the essential spectrum encroaches the right 
half A-plane. There are three types of instability [16 
3H3S]. The instability is called transient (TI) if the 
essential spectrum encroaches the right half plane and ab- 
solute spectrum are either in the left half plane or does 
not exist. The propagating DW emits stern waves shown 
in Fig. [3jb)i. The instability is called convective if both 
essential and absolute spectrum encroaches the right half 
A-plane. In this case, the emitted waves can propagate 
in both direction as shown by Fig. [3][b)ii. For an convec- 
tive instability, if any branching point is also in the right 
half A-plane, the instability is called absolute. An abso- 
lute instability can then emit non-traveling (zero group 
velocity) waves as illustrated in Fig. [3jb)iii. For LLG 
equation, since the absolute spectrum is the segment con- 
necting two branching points Sdi and Sd 2 (Fig. [3] (a)), 
the absolute instability (AI) and convective instability 
(CI) co-exist. It is known that transient instability is 
very weak that can be removed under proper mathemat- 
ical treatment [16|, [26j. Thus, we should not expect to 



have great physical consequences. On the other hand, 
the absolute instability move with the DW, and cause 
the change of DW profile [IiH2(| . It is known @ that 
field-induced DW propagating speed is proportional to 
the energy damping rate that is sensitive to DW profile. 
Therefore absolute instability, which deform propagat- 
ing DW profile, shall substantially alter DW speed. This 
may explain why the field-induced DW speed start to de- 
viate from the Walker result only when the field is large 
enough to emit both stern and bow waves in simulations 
10]. 

Fig. (3Jc) is the calculated phase diagram in K± and 
p = H/H c plane. A transition from transient instabil- 
ity (denoted as TI in the figure) to absolute/convective 
instability (AI/CI) occur at a critical field H% as lng as 
K± > K*[ w 0.085 at which H 2 = H c . It means no 
absolute/convective instability exist for K± < K^, and 
one shall not see noticeable change in famous Walker 
propagation speed mentioned early. This may explain 
why many previous numerical simulations on permal- 
loy, which have small transverse magnetic anisotropy, are 
consistent with Walker formula. A snapshot of the con- 
verting wavepackets could be identified in Fig. 2 in Ref- 
erence jlOj where wavepackets can be seen in the vicinity 
of the traveling DW and travel to both directions. 

It should be noticed that the effects of point spectrum 
have not been analyzed. In principle, it can also affect the 
stability of the Walker solution, and should be a very in- 
teresting subject too. Unfortunately, there are not many 
theorems on the point spectrum yet. Thus, one can only 
rely on a numerical method to find a point spectrum of 
operator L and to find out whether it can also induce any 
instability on a propagating DW. 

In conclusion, we showed that a Walker propagating 
DW will always emit stern waves in a low field, and both 
stern and bow waves in a higher field. Thus the exact 
Walker solution of LLG equation is not stable. The true 
propagating DW is always dressed with spin waves. The 
emitted spin waves shall be damped away during their 
propagation, and make them hard to be detected in real- 
istic wires. For a realistic wire with its transverse mag- 
netic anisotropy larger than a critical value and when 
the applied external field is larger than certain value, 
a propagating DW may undergo simultaneous convec- 
tive and absolute instabilities. As a consequence, the 
propagating DW will not only emit both spin waves and 
spin wavepackets, but also change significantly its profile. 
Thus, the corresponding Walker DW propagating speed 
will deviate from its predicted value, agreeing very well 
with recent simulations. 

This work is supported by Hong Kong RGC Grants 
(604109 and RPC11SC05). 
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